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Abstract. Wilson chains, based on a logarithmic discretization of a continuous spectrum, are widely used to 
model an electronic (or bosonic) bath for Kondo spins and other quantum impurities within the numerical 
renormalization group method and other numerical approaches. In this short note we point out that Wilson 
chains can not serve as thermal reservoirs as their temperature changes by a number of order AE when a 
finite amount of energy AE is added. This proves that for a large class of non-equilibrium problems they 
cannot be used to predict the long-time behavior. 

PACS. 71.10.-w Theories and models of many-electron systems 
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In 1975 Wilson pioneered [T] the "numerical renormal- 
ization group" (NRG) approach to solve the Kondo model. 
This approach is able to provide numerically exact solu- 
tions to the Kondo problem, probably the most important 
and most influential quantum impurity model. The NRG 
is by now widely used for a large range of fermionic and 
bosonic quantum impurity problems both at zero and fi- 
nite temperature 2 . 

A quantum impurity problem is defined by some local 
degree of freedom coupled to a (most often non-interacting) 
bath. In the case of the classical Kondo model this bath 
is given by an infinite system of electrons with a constant 
density of states at the Fermi energy. 

Wilson's numerical solution consists of two main steps. 
First the bath is discretized using energies defined on a 
logarithmic grid. 
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This logarithmic grid is optimized for the logarithmic renor- 
malization group flow characteristic for the Kondo model 
and allows to capture even the exponentially small energy 
scales characteristic for the Kondo problem. In praxis [2 
it is implemented via the so-called Wilson chain, a chain 
of local sites coupled by nearest neighbor hopping rates 
which decay exponentially in the distance from the quan- 
tum impurity. The second step, the actual renormalization 
group transformation, strongly builds on this logarithmic 
discretization but it plays no role for the following dis- 
cussion. 

An exciting recent development is that it is now pos- 
sible to study quantum impurity problems numerically 
in non-equilibrium using NRG [SHHOinilZ] • Especially F. 
B. Anders and A. Schiller pioneered powerful methods to 
study quantum quenches [1] and to investigate steady- 
state transport using many-particle scattering states ,5 . 



Also other numerical methods, most notably time-dependent 
density matrix renormalization group have been applied 
to Wilson chains "W . 

NRG and the use of Wilson chains has been extremely 
successful [2]. Therefore it is a relevant question to un- 
derstand whether there are limitations to this approach. 
In this short note we briefly discuss a very severe limita- 
tion: Wilson chains can not be used as thermal reservoirs. 
We argue that this implies that for certain classes of non- 
equilibrium problems Wilson chains give wrong results. 

A reservoir is by definition an infinite system which can 
be used to fix the chemical potential and the temperature. 
More precisely, when adding a finite amount of particles 
or a finite amount of energy, the chemical potential and 
the temperature do not change. While in experimental 
systems reservoirs are never really infinite, in praxis the 
large number of electrons, for example in a metallic con- 
tact attached to a quantum dot, provides an almost ideal 
realization of a reservoir. 

Within this definition, the Wilson chain can serve as a 
reservoir for particles as one can check by calculating ^ , 
for the non-interacting Wilson chain in the absence of the 
quantum impurity 
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Here N is the total number of particles, /(e) = 1/(6^' '^)/^-| 
1) is the Fermi function and the factor 2 accounts for the 
spin degeneracy. ^ diverges with the number of sites in 
the Wilson chain implying that in practical calculations 
the chemical potential stays constant when, for example, 
N increases by 1. 
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In contrast, the total energy of the Wilson chain and 
its derivative vifith respect to temperature T remain finite 
in the limit A^sites — ^ oo for finite discretization parameter 
A>1. 

dT- ^Z.T^^^"^^V-oo (l + e^-")2 - In^ 

(3) 

In practical calculations A is never very close to 1 (an 
often used value is yl = 2) therefore dE/dT is a number 
of order l! 

This has severe practical implications for non-equili- 
brium situations. Consider, for example, a quantum quench 
where at time i = suddenly a hybridization of strength 
r is switched on, which couples an Anderson impurity to 
a Wilson chain at T = 0. The following time evolution 
conserves energy. The ground state energy of the system 
with finite hybridization is approximately lower by F than 
the energy of the initial state. Assuming that the system 
equilibrates, Eq. ([3]) implies that final state has a huge 
temperature of order _rinyl, much larger than the Kondo 
temperature. In contrast, an Anderson impurity coupled 
to a real reservoir, where dE/dT diverges with system size, 
relaxes to the T = groundstate in the thermodynamic 
limit. While in real systems the energy can be transported 
away from the quantum impurity to infinity this is not pos- 
sible for the Wilson chain due to its finite heat capacity. 
Our argument proves that in the long time limit a wrong 
result is obtained for all quantum quenches where the en- 
ergy of the initial configuration is higher than the ground 
state energy by an amount of the order of the energy scale 
one is interested in. With Eqn. Q one can estimate quan- 
titatively the heating of the quantum impurity. This can 
serve as a quantitative estimate of the backreaction to the 
quantum impurity which arises as the Wilson chain is not 
able to transport energy away to infinity. 

It is less clear to what extent the short-time dynamics 
is affected but even in this case it would not be surpris- 
ing if the problem to transport energy away is of relevance 
even shortly after a quantum quench. Recently, P. Schmit- 
teckert [S] made a closely related observation by studying 
the time-evolution of a wave packed on a Wilson chain: 
the wave packed piles up in a "Wilson Tsunami" and is 
not able to reach the end of the chain. 

It is an interesting question whether and to what ex- 
tent the finite heat capacity of Wilson chains is also of im- 
portance when using scattering states to treat steady-state 
non-equilibrium [5]. Consider, for example, a quantum dot 
in the Kondo regime in the presence of a finite bias voltage 
[TD] . In this case the finite product of voltage and current 
implies that power is pumped into the system. This power 
is ultimately dissipated in infinite reservoirs where the pre- 
cise form of the dissipation mechanism is usually not im- 
portant as the relaxation occurs far way form the impurity. 
In the case of an infinite Wilson chain where energy can- 
not be added without modifying the bath, one can suspect 
that (as in the case of a quantum quench) the interacting 
non-equilibrium situation cannot be treated properly. In- 
terestingly, present NRG implementations seem not to be 



able to recover the result pUfTU] that the Kondo peak in 
the spectral function splits in the presence of a bias volt- 
age which exceeds the Kondo temperature. This result 
is well established by (renormalized) perturbation theory 
and controlled one-loop renormalization group [T^ calcu- 
lations. Possibly, the numerical results which have been 
obtained for the Anderson model where not sufficiently 
deep in the Kondo regime to see the splitting. Alterna- 
tively, it is also possible that the use of Wilson chains is the 
origin of the problem. Note, however, that also other nu- 
merical approaches which are not relying on Wilson chains 
have presently problems to obtain the splitting [T^ . 

In contrast to the fermionic Wilson chain, for bosons 
the heat capacity of the Wilson chain is infinitely large 
due to the possiblity to put a large number of bosons into 
the states with very low energy (replacing the plus by a 
minus sign in the integral of Eq. ([3]) leads to a divergence). 
In practical NRG implementations [2], however, only a fi- 
nite number of bosons is taken into account which renders 
dE/dT again finite. 

In conclusion, it seems that the use of Wilson chains for 
non-equilibrium problems is a much more dangerous ap- 
proximation compared to the equilibrium case as Wilson 
chains cannot take up finite amounts of energy without 
changing their properties. 

I thank F. B. Anders, R. Bulla, Ch. Schiitte, E. Sela, K. Ro- 
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608 of the DFG. 
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